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Abstract 

We calculate the nucleon and the delta mass to fourth order in a covariant formulation of the 
small scale expansion. We analyze lattice data from the MILC collaboration and demonstrate that 
the available lattice data combined with our knowledge of the physical values for the nucleon and 
delta masses lead to consistent chiral extrapolation functions for both observables up to fairly large 
pion masses. This holds in particular for very recent data on the delta mass from the QCDSF 
collaboration. The resulting pion-nucleon sigma term is a-j^N = 48.9 MeV. This first quantitative 
analysis of the quark-mass dependence of the structure of the A(1232) in full QCD within chiral 
effective field theory suggests that (the real part of) the nucleon-delta mass-splitting in the chiral 
limit, Aq = 0.33 GeV, is slightly larger than at the physical point. Further analysis of simultaneous 
fits to nucleon and delta lattice data are needed for a precision determination of the properties of 
the hrst excited state of the nucleon. 
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1 . The A(1232) is the most important baryon resonance. It is almost degenerate in mass with the 
nucleon and couples strongly to pions, nucleons and photons. It was therefore argued early that 
spin-3/2 (decuplet) states should be included in baryon chiral perturbation theory which is the 
low-energy effective field theory of the Standard Model. In Ref. ^ and subsequent works use was made 
of the heavy baryon approach, which treats the baryons as static sources. However, due to the fact 
that the nucleon-delta mass splitting A = ttia — misf stays finite in the chiral limit and thus introduces 
an additional low energy scale, in chiral effective field theories (ChEFT) special care has to be taken 
about the decoupling of resonances in the chiral limit [ 2 ]. This requirement can be systematized by 
counting the nucleon-delta mass splitting as an additional small parameter [2], which ensures that 
at each order in the chiral expansion enough counterterms are present to guarantee decoupling and 
renormalization |^. The corresponding power counting was called the “small scale expansion” (SSE) 
|2]. The heavy baryon approach has been successfully applied to a variety of processes, for reviews see 
and a status report on chiral effective field theories with deltas is given in [2j. More recently, it 
was realized that for certain considerations/processes a Lorentz-invariant formulation of baryon chiral 
perturbation theory is advantageous. A particularly elegant scheme to perform covariant calculations 
is the so-called “infrared regularization” (IR) of [S]. The renormalization of relativistic baryon chiral 
perturbation theory has been discussed in detail in [Hj. In Ref. m we gave a consistent extension of 
the infrared regularization method in the presence of spin-3/2 fields. It was in particular shown that 
in the covariant formulation of the SSE the contribution of the non-propagating spin-1/2 components 
of the Rarita-Schwinger field can be completely absorbed in the polynomial terms stemming from the 
most general effective chiral Lagrangian, analogous to the situation in the non-relativistic SSE [3]. In 
this letter, we apply the covariant SSE formalism to the nucleon and the delta mass"^® as well as the 
ttN sigma term. We perform a fourth order calculation in the small parameter e, where £ collects 
small external momenta, the pion mass and the NA mass splitting. These explicit representations of 
m-AT and ttt-a serve as chiral extrapolation functions to analyze lattice simulations for nucleon and delta 
masses involving dynamical fermions^®, as for example reported by the CPPACS |14] and JLQCD 
collaboration Hg. While a detailed (combined) analysis of these and forthcoming data jib) analogous 
to the work reported in Ref. m certainly also needs to take into account the observed finite volume 
dependences nulls], in this letter^^ our aim is more modest. We attempt to connect the lattice data 
for nucleon and delta masses as reported by the MILC collaboration nsj with the known results at the 
physical point. The MILC data have the exciting feature that they cover unphysical pion masses as 
low as — 350 MeV on a relatively large lattice with L = 2.6 fm, although there is still considerable 
discussion within the lattice QCD community about the technical details of the employed (staggered) 
fermion actions, see e.g. noj. Clearly, such extrapolation functions based on chiral perturbation theory 
cease to make sense at too large values of the quark (pion) masses, but as we will demonstrate later, 
we can nicely capture the trend of these data^®. More lattice results also at lower pion masses are 
clearly needed, and these should be analyzed utilizing the (infinite volume) extrapolation functions 

’^®For the substantial literature on this topic analysed within the heavy baryon approach we refer to the most recent 
work CD and the literature cited therein. For related early work on the problem of chiral extrapolation functions for 
A (1232) properties see 1121 . 

’^^For a discussion of the chiral extrapolation of nucleon and delta masses in quenched QCD see m 

’^’^Some of our results have been reported previously m 

’^^The bulk of the MILC data have been obtained in simulation runs with three active flavours, while our chiral 
extrapolation functions apply to a scenario with only two light quark flavours. However, at the (still relatively large) 
quark masses studied in Ref. ca. corresponding to pion masses > 350 MeV, the MILC collaboration reported no 
significant differences between their two and three flavour runs for the observables considered here. Note that they 
published a single two-flavor run at one quark mass that showed no sizable deviations from the three-flavor runs in m 
We will therefore treat the MILC data as if they constituted dynamical two flavor results. This is further corroborated 
by our analyses of the nucleon mass in SU(2) and SU(3) heavy baryon CHPT - no large effects were found when going 
from two to three flavors, see m and m respectively. 
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given here. 


2 . Our calculations are based on the effective Lagrangian of nucleons and deltas coupled to external 
sources. The various contributions to S-matrix elements and transition currents are organized in 
powers of the small parameter e, where e collectively denotes small pion four-momenta, the pion 
mass, baryon three momenta and the nucleon-delta splitting, A = (more precisely, the 

difference in the chiral limit). The expansion in powers of £ is called the small scale expansion. In 
what follows, we will consider the nucleon and the delta mass to fourth order in the SSE. Consider 
first the pion-nucleon Lagrangian. The terms pertinent to the observables calculated here read (for 
details, see 001221 ), 
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Here, i/’A'' denotes the nucleon spinor, is the axial-vector coupling constant (in the chiral limit), 
mo the chiral limit value of the nucleon mass, x+ contains the external scalar fields which contains 
the quark mass matrix and (...) denotes the trace in flavor space. The Cj are low-energy constants 
(LECs). Their values can e.g. be determined in the analysis of pion-nucleon scattering in the delta-full 
EFT |24| . It is known that in particular C 2 and C 3 are largely saturated by A-exchange [25| . thus the 
remaining contributions in an EFT with explicit spin-3/2 fields are expected to be very small. The 
counterterms B 23 , B 32 , Ei, E 2 in appear naturally in SSE (Aq ~ 0{e)) and are required for 

the renormalization of the loop graphs at 0(e^’p. Furthermore, their finite parts are chosen in such 
a way that in the limit /S.o 0 the 0{e'^) SSE result recovers the 0{p^) result of heavy baryon 
chiral perturbation theory (HBCHPT). The fourth-order LECs Cj contribute to the nucleon mass in 
the combination ei = Ibeas + 2 eii 5 + 6110/2 (for details, see [IHII^). The relevant terms from the 
covariant ttNA Lagrangian read m 


= CAi^aO^A w'g'ipN + h.c. , (5) 

+ —+ + , ( 6 ) 

I "2-0 J 

with re/ = (r*tt^)/2 and = (r* Ui/])/2, where r*, i = 1,2,3 denote the Pauli matrices in 

isospin space. The tensor = gag + takes care of the “point-invariance” of the spin- 

3/2 theory jS]. Here, "i// denotes the spin-3/2 field in Rarita-Schwinger notation with an additional 
isospin-index i and the constraint t* "f// = 0 [ 2 ] . The LEG ca represents the leading axial-vector N A 
coupling constant (frequently called g-jrNA in the literature). At second order, we have two relevant 
couplings, parameterized by the LECs 63 and bo. We note that in the non-relativistic theory the 
leading contribution of bo can be absorbed in the coupling 63 to O(e^) |21. In the non-relativistic 
SSE bo therefore only appears as an independent coupling at O(e^) . Finally, we list the terms of the 
Lagrangian that describe the delta propagation and its coupling to pions and external sources. The 
pertinent covariant structures up to 0 (£^) have been constructed in 1221 , here we add further terms 
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needed up O(e^) in close analogy to the nucleon case (see also HH). Specifically, we employ 
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with {t^D^), and denotes the isospin 3/2 projector, ^*7 = 

jd _ rVAS- Further, gi denotes the axial A A coupling constant (in the chiral limit), whereas m\ 
refers to the A(1232) mass (in the chiral limit). The LECs a, correspond^® to the c* in Eq. © and so 
on. The Bf" and Ef^ are taken such that the values of tuq and ai in the chiral limit are not affected 
by loop corrections. 


3 . Choosing A = —1 the propagator for a spin-3/2-isospin-3/2 particle in the Rarita-Schwinger 
formalism has the general form"^^® 
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The Dirac-tensor S^j^y can be written as a linear superposition of spin-3/2 and spin-1/2 projection 
operators P%^ 
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For details regarding the properties of the projectors we refer to m- As can be clearly seen from 
Ea. dT^ only the spin-3/2 components are associated with a propagation, whereas the (spurious) spin- 
1/2 components correspond to local contact operators. Given that the chiral effective field theory 
for a coupled nucleon-delta system contains the most general set of local contact operators allowed 
by chiral symmetry, it was argued in cni that in chiral EFT one only needs to take into account 
the propagation of the spin-3/2 degrees of freedom. The effects of the spurious/off-shell spin-1/2 

note that while only three of the d LECs are contributing to the mass functions considered here, there are 
actually 9 different LECs ai that should be taken into account. However, it can be shown that there are only three 
different linear combinations of them contributing to the masses. We therefore take ai, 02 , 03 as representatives of these 
three independent structures. 

’^'^°The small scale expansion constitutes one possible form of an effective low energy theory of QCD for low lying baryon 
resonances. Its range of validity and applicability therefore clearly lies in the low energy domain with momenta below 
the chiral symmetry breaking scale A,,.. All physics connected with higher energies, e.g. the field theoretical deficiencies 
of the Rarita-Schwinger approach discussed in the literature m, is therefore mandated to be accounted for only via the 
counterterms/local operators of the theory. Gauge invariance can be shown to hold perturbatively for all forms of the 
propagators addressed in this work, see Ref. 1281 . For alternative approaches to spin-3/2 particles in ChEFT see 12911301 . 
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components are already completely accounted for by the counterterms of the theory. In m we have 
therefore proposed to utilize the propagator form A 


with the spin-3/2 projector 
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In the meantime it was realized in practical applications that it is advantageous to use a different form 
for the propagator than the one given in En. (j13jl (see also the related work in [31 jl. In Ref. |1 D] it was 
shown that due to the extra 1/p^ structure in Ea. (ll4jl the standard rules of infrared regularization had 
to be modified for a consistent treatment of the low-energy modes, making loop-calculations in this 
covariant field theory more tedious. Here we propose another separation of the general Dirac-tensor: 
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is not connected with the spin-3/2 propagation, that it is local and that it cannot generate any 
contributions which are not already contained in the counterterms of the effective field theory. As an 
alternative to Eq. m we therefore propose"^^^ the propagator form 
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We note that in form B there is no more 1/p^ structure, making the evaluation of loop diagrams 
containing spin-3/2 degrees of freedom in the IR-scheme a lot simpler. However, we emphasize that 
in a low-energy (chiral) effective field theory the finite/renormalized results obtained via En. (n3l) or 
Eg. ([1711 are identical"^^^—it is merely a matter of practicability which method is chosen for a particular 
observable. In what follows, we will work with form B. We note explicitly that one could also work 
with the full Rarita-Schwinger propagator of Eq. (HU), however, form A or form B have the advantage 
that the results can be directly mapped to the known couplings of the heavy-baryon formulation of 
SSE where also only the spin-3/2 components are part of the propagator ^IB]- Finally we remark that 
for low energies the width of A(1232) can be treated perturbatively in the covariant SSE, analogous 
to the situation in the non-relativistic theory. Throughout this work we will therefore not address 
this issue, applications of the covariant SSE in the resonance region will be discussed in a forthcoming 
paper m 

4. Armed with the effective Lagrangian, we are now in the position to calculate the nucleon and the 
delta self-energy to fourth order in the covariant SSE. The covariant formalism allows in particular to 
resum the contributions from kinetic energy insertions. The nucleon mass can be written as 


m = mo - 4 ciM2 -4M2AoH 23 - HsaA;] -4ei(A)M^ - EiA^ - dM^A^Ea +mN-^°°P + m^-^°°P , (18) 

note that this form of the spin-3/2 propagator is also employed when the interactions of the spin-3/2 field are 
constructed according to the rules of gauge symmetry for a massless spin-3/2 field, e.g. see HO]. 

Obviously the high-energy behavior between the two propagator expressions is different, leading to different numerical 
values for the associated counterterms to yield the same total result. 
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Figure 1: Fourth order contributions to the nucleon (left panel) and the delta mass (right panel). Solid, 
double and dashed lines denote nucleons, deltas and pions, in order. The corresponding dimension 
two {ai,bi,Ci) and four (e*) LECs are also given at the pertinent vertices. 


with the corresponding fourth-order tree and one-loop graphs shown in Fig. ^ The contributions from 
the nucleon and delta loop graphs, which start at 0{e^), take the form 
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Here, is the pion decay constant and A the scale of dimensional regularization. Throughout, we 
set A = 1.232 GeV. The loop functions which appear here and for the delta mass discussed below are 
defined as: 
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with 
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Using these, the A-loop contribution can be written as an infinite series in (modulo chiral logs): 

m^-'°°P = /o + /iM2 + /2M4 + ... , (29) 

with lo, li being infinite series in Aq starting at O(e^) and 0{e), respectively. Note that the decoupling 
theorem is automatically fulfilled here since there is no nonanalytic contribution ~ M^. Now as a 

minimal way of determining the LECs B 23 , i? 32 j Ei and E 2 one can choose B 23 + AqE 2 = /i/4, and 

i ?32 + AgEi = lo/4:, which means that the A contribution will first occur at order M^. 

The chiral expansion of the delta mass can be written in a way similar to the case of the nucleon 
Eq. (jlSf) (the corresponding tree and loop graphs that contribute at 0{e^) are collected in Fig. P), 

rriA = m%-4aiM^-4M^AoB^-B^Al-4et{X)M^-E^A^-4M^AlE^+ml-^°°^+m^-^°°P , (30) 
with 

^2 


32tt'^P^ V ^ _ 77j2 _|_ ]\^2 _ yj\(p^^rri^^ M‘^) 

if X{p‘^, m'^, M^) > 0. Otherwise, we have 


Jo — 


1 


—X{p^,m?^M^) arccos — 


P 


^ — m? Ml ' 


167r2p2 

The other relevant loop functions can be expressed as 
r(bl) 7.„2 _2 


2 M^V^ 


Tl 


(/, ^ (^(p 2 -rr? + 




4(J — l)p2 


A(p^, m^, Ml)lff^ - {p^ -m^ + Ml)A, 


Tl 


= ^((P^ -m^ + Ml)l[lJ'> + 


All 


We also employ the standard definition of Ajr, 


d'^k 


1 


m 


N—loop C 

F 2 

+ bs 


caMI 


+ 


^ {(m2, + mo) Jn°^((m2,)^m^) - ^^((m2,)^mg) } 

|a^ + [(m2,)^ -ml- Ml]l[l'°\{m%)^,ml) + 2mom2,/{l’^^((m2,)^, mg)| 
|A^(mo + m2,) + [(m2,)^ - mg + Ml] 


FI d 
be CA Ml 
naFl d 


(mo + m2,)/^2’°^((m2,)^mg) - m2,/lg’"^((m2,)^ mg) 


0 \2 2^ 


(31) 


and 


m 


A—loop _ 


bhl 


6 (m 2 ,) 2 F 2 

1 


m 2 , 




^ -M 2 (m 2 ,)^/{?’°^((m 2 ,)^(m 2 ,)^) 


+ d^(4("^A)^ + {d- mmrum^Ar, (^ 2 ,)^) 


2 lr( 2 - 0 )/. 0 ^2 / 0 n 2 ^ 


+ m 2 . 




J — 1 


Ml 


A, 


1 


d{d-l)^"^A El 


3^' ^_2ai + :^ + B) A^ , 


(32) 


7 





























Figure 2: The nucleon mass (dot-dashed line) and the (real part of the) delta mass (solid line) as a 
function of the pion mass. The filled diamonds denote their physical values at the physical pion mass. 
The dashed line is the chiral extrapolation for the A based on SU( 6 ) as explained in the text. The 
filled squares and circles are the MILC data m- The filled triangles are the recent data from QCDSF 

m- 


with A = a 2 + ■■ ■ and B = + .... It can easily be verified that starts as whereas, as 

in the case of the nucleon mass, has a constant and a piece Ml 

^N-loop ^ ^ ^ ^ . ^ 33 ) 

As before, we will choose B^, in such a way that starts to contribute only at 

Ml that is B^ A^E^ = li/4 and B^ + AqE^ = Iq/4. Note that we did not take into account the 
difference between the chiral limit value of the pion mass and its physical value m- This effect is 
known to be small, see e.g. the three-flavor calculation [ 22 j . 

5. We are now in the position to analyze the nucleon and delta mass formulae given in Eqs. (j1 SIHOjl . 
They contain a certain number of LECs, some of which are (not very accurately) known from the study 
of pion-nucleon scattering in the heavy baryon SSE m Here our aim is modest^^^: In addition to 
the known values at the physical point we take the data from MILC m for the nucleon and the delta 
as function of the pion mass and try to describe these with LECs of natural size. Such a description 
is indeed possible, as shown in Eig. So we do not intend detailed least-square hts here but rather 
try to find out whether the existing data shown in the this figure can be consistently described by 
our mass formulas with LECs of natural size. We stress again that a more refined analysis of e.g. 
pion-nucleon scattering in the covariant SSE is mandatory to put stringent constraints on certain 
combinations of the LECs (see also the extensive discussion in Refs. ED ESI on this issue). We have 11 
(combinations of) parameters to determine, these are ci, C 2 / 4 -I-C 3 , ai, ei, ef^, ca, hA, 63 , be, A/4+B and 
Aq. (i) constrained parameters: ca and Aq are hxed from the mass and width of the delta. We note 
that the NA mass splitting in the chiral limit, Aq = 0.33 GeV, indicates a slightly larger NA mass 
splitting in the chiral limit than at the physical point. A similar observation was made in the case of 

Ultimately global, simultaneous fits of several nucleon and delta observables calculated within covariant SSE to 
next-to-leading one-loop order need to be undertaken to obtain reliable information on LECs from lattice QCD. 
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quenched QCD, see m (for discussion about the difference between the mass splitting in the quenched 
approximation and full QCD, see also IMI). The LECs from the pion-nucleon Lagrangian, Eqs. (EH 
can be estimated from pion-nucleon scattering in the presence of the delta. We use ci = —0.8GeV“^ 
(which is within the uncertainty of the values determined in e.g. m) and Cl = C 2 = 0, C 3 = 0.5 GeV ^. 
The small values of € 2,3 are consistent with resonance saturation studies of EHI and the fits in M- (ii) 
Fit parameters: The remaining parameters are fit to the masses. We find for the two axial dimension 
two NA LEGs the values 63 = 0.75 GeV“^ and = —0.75. The axial AA coupling is found to be 
Ha = 2, which is not far from the SU( 6 ) or large-Ac value Ha = 9gA/5 = 2.28. Furthermore we get 
ef' = —lGeV“^, ai = —0.3GeV“^, A = 0 and B = 0.5GeV“^. These are all natural values. It is 
interesting to note that ai is markedly smaller than ci, although both couplings should be equal in 
the SU( 6 ) limit. We refrain here from a detailed study of the theoretical errors that will be given in 
a forthcoming publication. Still, it is interesting to study the strict SU( 6 ) limit. In that case, one 
would have ai = ci = —0.8 GeV“^ and Ha = 2.28. As can be seen from the dashed line in Fig. [21 the 
assumption of strict SU( 6 ) symmetry is clearly at odds with the MILC data, indicating that ai and ci 
indeed seem to have different values. Also shown in Fig. [21 are the recent QGDSF data for ruA, which 
were not used in the fit but are nicely consistent with our extrapolation function. Note also that the 
QGDSF data are based on two-flavor simulations and are not very different from the MILG data in 
the region of overlap. This further supports our assumption on the treatment of the MILC data. We 
stress again that the resulting values of the LECs are to be considered indicative and a more analysis 
employing also contraints from other physical processes should follow. 

From the small value of the LEC oi one immediately deduces that the ttA sigma term appears to 
be significantly smaller than its nucleon cousin because at leading order in the quark mass expansion 
we have + ... and ct^a = —-|- .... It is clear that this interesting observation 

deserves further study. Finally we note that the sigma term for the nucleon resulting from this “rough” 
fit is found as 

= 48.9 MeV , (34) 

to order e^. We note that this value is consistent with the classical result of Ref. m, which was 
confirmed in inn] in a HBCHPT analysis of pion-nucleon scattering and in m in a ChEFT analysis 
of lattice data in a formalism without explicit delta degrees of freedom. It is also in agreement with 
the recent CHPT analysis of the three-flavor MILC data, see m- For the vrA sigma term we get 
o^ttA = 20.6 MeV. Again, these results need to be refined and bolstered by more detailed precise fits 
to the lattice data including also error and correlation analysis including also lattice data on other 
observables - the mass data alone are not sufficient to precisely pin down all parameters. Such an 
analysis, however, goes beyond the scope of this paper. 

6 . In this letter, we have presented a covariant extension of the small scale expansion, extending 
our earlier work m- We have analyzed the nucleon and the delta mass in view of the lattice data 
from MILC. These data seem to indicate sizeable SU( 6 ) breaking effects that can be interpreted as 
a reduced pion cloud contribution in the resonance field (the delta) as compared to its ground-state 
(the nucleon). Our results obtained from fits to full QCD lattice data further indicate that the NA 
mass splitting can be larger in the chiral limit than at the physical point, analogous to observations in 
quenched QCD. We have also calculated the pion-nucleon sigma term and found good agreement with 
earlier determination from pion-nucleon scattering and the analysis of SU(2) and SU(3) lattice data 
for the nucleon. The recent QCDSF data for ruA are also consistent with our chiral extrapolation. It 
would be very valuable to have more (2 flavor) lattice results for the delta preferably at small pion 
masses to further sharpen these conclusions. 
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